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Abstract. In molecular crystals like McPTCDI and PTCDA the molecule 
are regularly arranged creating quasi-one-dimcnsional molecular stacks. The 
intermolecular distance in a stack (about 3.3 A) is comparable with the electron- 
hole distance in the excited molecule. The mixing between Frenkel excitons (FEs) and 
charge-transfer excitons (CTEs) is very essential for the excitonic and vibronic spectra 
of both crystals. In this paper, we make simulations of the linear absorption spectra of 
the abovementioned crystals. The basic Hamiltonian describes the FE-CTEs mixing 
in the molecular stack (point group d) caused by two transfer mechanisms, notably 
of the electron and the hole on the neighbor molecules. The vibronic spectra consist 
of mixed excitons and one vibrational mode of an intramolecular vibration linearly 
coupled with FE and CTEs. Using the vibronic approach, we calculate the linear 
optical susceptibility in the excitonic and one-phonon vibronic regions of the molecular 
stack, as well as of a crystal which contains two types of noncquivalent stacks. We put 
the excitonic and vibrational parameters for the crystals of PTCDA and MePTCDI 
fitted in the previous studies. We analyze the general structure and some important 
features of the lineshape of the linear absorption spectra in the spectral region of 
15 000-23 000 cm -1 . We vary the values of the excitonic linewidth and the parameters 
of the linear exciton-phonon coupling and look for the values which reproduce the 
absorption lineshape similar to the absorption spectra of the investigated crystals. Our 
study exhibits the necessity of introducing the FE-CTEs mixing in the interpretation 
of the linear absorption spectra, especially of the MePTCDI crystal. 
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1. Introduction 

The mixing of Frenkel excitons (FEs) and charge-transfer excitons (CTEs) has been 
studied both experimentally and theoretically in many one-component molecular 
crystals, e.g., polyacenes [U El 131 E], perilene derivatives [SI El [7], fullerenes [8j [9], 
and other. The excitonic and vibronic spectra of quasi-one-dimensional crystals like 
3,4,9, 10-perylenetetracarboxylic dianhydride (PTCDA) and N, iV'-dimethylperylene- 
3,4,9, 10-perylenetetracarboximide (MePTCDI) have been treated in (5j EJ [7j and 
PHI [HI [121 EH E] • The short molecular distance in quasi-one-dimensional stacks causes 
a strong FE-CTEs mixing and mixing of their vibronic satellites as well, especially in 
the absorption spectra in the spectral region of 2-3 eV (15 000-23 000 cm -1 ). 

In the present paper, we calculate the linear absorption spectra of one-component 
MePTCDI and PTCDA crystals applying the vibronic approach developed in our 
previous papers [131 US]- It allows complex calculations of the absorption in pure 
excitonic, one-phonon vibronic, two-phonon vibronic, etc. spectra. In paper [16] the 
vibronic approach is the tool of studying the FE-CTEs mixing in a two-component stack 
of alternatively arranged donor-acceptor (DA)-molecules. While the FE-CTEs mixing 
in [16] is a probable but hypothetical model, in the present study we turn to the real 
absorption spectra and use the parameters of the FE-CTEs mixing from Refs. [51 l6l [TO] . 

Two general differences exhibit the models' comparing of the FE-CTEs mixing in 
one- and two-component molecular stacks [16], notably: (i) The inter molecular transfer 
both of electrons and holes must be a feature of the model for the case of one-component 
stack whereas only one transfer mechanism is sufficient in a DA-molecular stack. In the 
last case the Frenkel exciton represents a collectivized electronic excitation of a donor or 
an acceptor and the strongest transfer on the closest neighbor would be either of a hole 
or correspondingly of an electron. Obviously in an one-component stack both types of 
transfer on the neighbors could be of the same probability. The two-step processes of 
successive transfer of the electron and the hole on the same molecule ensure the transfer 
of a FE even in the case of its relatively weak direct intermolecular transfer [10J. In 
this way, our present study extends the calculations in [13] where only one transfer 
mechanism has been considered, (ii) In the most widely studied DA-crystal, anthracene- 
PMDA the excitonic absorption lines are narrow [T71 [THl [19] and many details in the 
vibronic spectra can be seen. In PTCDA and MePTCDI crystals the absorption lines 
are two order of magnitude wider. Thus we pay attention to the general structure of the 
excitonic and one-phonon vibronic spectra supposing a width of the excitonic lines of 
300-500 cm -1 (not 2-10 cm -1 which is the absorption width of the anthracene-PMDA). 

In the next section of the paper we involve the initial Hamiltonian in the case of 
a FE-CTEs mixing. The Hamiltonian contains one mode of intramolecular vibration 
linearly coupled to the FE and CTEs. In section 3 the linear optical susceptibility 
has been calculated in the excitonic and one-phonon vibronic spectra. In section 4 the 
linear absorption has been modelled using the excitonic and vibrational parameters of 
the PTCDA and MePTCDI crystals fitted in Refs. p [TO]. Section 5 contains some 
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conclusions. 

2. Hamiltonian for the case of a FE— CTEs phonon coupling 

We consider the excitonic and vibronic excitations in a linear molecular stack of TV 
identical molecules which are regularly arranged at a distance d each other. The 
point group of symmetry of the stack is Cj as is for the molecular stacks of PTCDA 
and MePTCDI [6J. The origin of the Frenkel exciton is a non-degenerate molecular 
electronic excitation with excitation energy Ep and L being the transfer integral between 
neighboring molecules. We denote by B n (B+) the annihilation (creation) operator of 
the electronic excitation on molecule n and get the following FE-part of the Hamiltonian: 

Hfe = ^2 EpB^Bn + 2^ L (S n ' <n+ i + 5„>-i) B+,B n . (1) 

n n,n' 

As usually, we consider two CTEs of equal excitation energy E c and C Uj i is the 
annihilation operator CTE, a — 1, with hole located on the site n and electron on 
the site n + 1, whereas the electron of the second CTE, a — 2, {C n ^) is located at 
molecule n — 1. We neglect the transfer of CTEs as a whole and the mutual coupling of 
the two CTEs since those processes can be realized through the transfer of the electron 
or hole at distance 2d which is less probable than the FE-CTEs mixing caused by the 
electron (hole) transfer at the neighbor molecule (see Refs. [El [IB])- We obtain the 
following CTEs-part of the Hamiltonian: 

H C t= E * C tC™ (2) 

n,o-=l,2 

and suppose the following operator for the FE-CTEs mixing: 

n 

+ e hl B+C n ^ 1A + e h2 B+C n+lt2 + he] , (3) 

where e e i and e e 2 are the transfer integrals of the electron from molecule n to molecules 
n + 1 and n — 1 correspondingly, and e^i, denote the transfer integrals of the hole 
from molecule n to molecules n + 1 and n — 1. Certainly the model with four transfer 
integrals is more complicated than the model in Refs. P, [7] with two mixing parameters 
only. But our model is more realistic because we take into account the inclination of 
the flat molecules of PTCDA and MePTCDI relatively to the stack axis. 

One intramolecular mode is only supposed to be coupled with the FE and CTEs 
and the phonon part of the Hamiltonian is 

-ffph = ^2 ^oa^n, (4) 

n 

where uq is the vibrational frequency and a n is the annihilation operator of one 
vibrational quantum on molecule n. The linear exciton-phonon coupling only is 
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manifested in the treated crystals [7] and we get the following exciton-phonon part 

^ex-phon = E [& B n B n + £C+ a C na ] hu (a+ + a n ) , (5) 

n,<r=l,2 

where £f and £ are dimensionless parameters of the linear FE-phonon and CTEs-phonon 
coupling, correspondingly. 

The full Hamiltonian H contains all the parts (CO)-© and can be transformed using 
the canonical transformation which eliminates the linear exciton-phonon coupling ([5]), 
see [201 El], 

H x = exp(Q)#exp(-Q), (6) 

where 

Q= E [tFB+B n + £C+,C nff ](a+-a n ). (7) 

n,a=l,2 

We introduce the vibronic operators 

V n = exp(Q)B n exp(-Q) 1 (8) 

U nta = exp(Q)C na exp(-Q) (9) 

and get the following transformed Hamiltonian 

#1 = E (^F - ^o£ F ) V+V n + Y, L (<*n',»+l + <Wl) V+V n 
n nn' 

+ J2( E c~ U+a U na + E ^a>n 

n,cr n 

+ J2 [teiV+Um + e o2 V+U n2 + e hl V+U n - 1A + e h2 V+U n+1 , 2 + h.c] (10) 

n 

In a stack with inversion center, point group Cj, the excitons are gerade or ungerade 
(in the center of the Brillouin zone, at k — 0). The ungarade excitons only are dipole- 
active and influence the linear optical susceptibility and the absorption spectra. In the 
case under consideration, the operator of the transition dipole moment has the following 
form pElE]: 

P = E [PF W + V ») + Pct (^2 - U nl + [/+ - [/+)] • (11) 

n 

The gerade FEs, as well as the symmetrical combination of CTEs, (U^ 2 + ^Ci) 
can be also mixed, but due to their vanishing transition dipole moment they will not be 
considered here. 

Introducing the Fourier transform in the momentum space of the vibronic operators, 

Vk = -7= ^ V n exp(ifcna), (12) 
v N — 

U k ,a = —j= E U Hta exp(ikna), (13) 

* n 
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we obtain the following form of the Hamiltonian fllOp 

k k,a 

+ Yl i ^ + £ 'h COS k + K Sin k \ V k ( U k,l + U k , 2 ) 
k 

+ [e" + e'i cos k + i^ sin k) (U k)2 - U ktl ) + h.c.} + ^ fku a+a n , (14) 

where 

< = (e e i + e c2 ) /2, e' h = (£hi + e h2 ) /2, (15) 
<' = (£ e2 - e e i) /2, e£ = (^ 2 - £hi) /2. (16) 

For the case of ungerade FEs, their mixing with the symmetrical (even) combination 
(£4,i + U k>2 ) a t k = is impossible, and thus the mixing parameters e' e and e:^ vanish, 
e' c = e' h = 0. The final expressions for operators (j3J) and ffTTl) are 

^fc = £ K< + £ h cos fc) (U kj2 - f/ M ) 

k 

+ ie'i sin fc^+ (C/ fe , 2 + U ktl ) + h.c] (17) 

and 



P = VN[pf (V k=0 + V k t ) 

+ PCT (U k=0 , 2 ~ U k=0jl + E/+ =0 ,2 - ^ + =0,l)] 



respectively. 



3. Calculation of the linear optical susceptibility 

The linear optical susceptibility can be calculated by using the formula [22] 

= H ^ u + ie ) + + ie )]} ( 19 ) 

with 

= -^(^(01^(0^(0) + P,-(*)Pi(0)|0), (20) 

where 1/ is the crystal's volume [in our case being proportional to Nv (v is the volume 
occupied by one molecule)] and P is the operator ffT8]) . The Green functions fl20l have 
been calculated as average over the ground state |0) only by taking into account the 
large values of E F , E c , kuo ^> k B T. 

We calculate the Green functions (I20p following the vibronic approach [TJ)J [T6] . In 
the next expression, the x axis is supposed to be oriented along the vector pp which 
includes angle 7 with the vector Pct ; and a = \pct/pf\- Then we can represent the 
linear optical susceptibility of one stack as 

p 2 1 99 

Xxx = — tt^t \ a 2 + 4aai2 cos 7 + 2a 2 «i cos 2 7] . (21) 

v a±a 2 — 2a( 2 L J 
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The functions «i, a 2} ctu have been calculated for the excitonic and one-phonon vibronic 
regions (see below). The PTCDA and MePTCDI crystals contain two types A and B of 
parallel molecular stacks, however, the excitonic and vibronic excitations in each stack 
interact very weakly with the excitations of the other stacks. In the same way as in 
Ref. [13], we calculate the crystal's susceptibility in an oriented gas model. We denote 
by 2tp the angle between the vectors pp and pf of two different stacks and suppose that 
these vectors are positioned in the (XY) plane, the crystal X axis been oriented along 
the sum pp + p F [23] . 

The components of the linear optical susceptibility of the crystal correspondingly 

are: 

PF ^ r 2 

Xxx = 7TT «2 cos ip 

+ 4a«i2 cos 2 Lf cos 7 + a 2 ct\ (1 + cos 27 cos tp)] (22) 

and 

Pf 2 r • 2 

Xyy = "2 sin p 

v ctict2 — 2af 2 

+ Aaa\2 sin 2 ip COS7 + a 2 «i(l — cos 27 cos </>)] . (23) 
We find the following expressions for functions a>i, a 2 , ol\ 2 : 

(1) In the excitonic region expressions practically coincide with formulas in Ref. |16j . 
namely 

ai = hw - (E F + 2L) - m 0F (l), (24) 

a 12 = e" c + el (25) 

a 2 = h[uj-VL Gc {l)} -E c , (26) 

where Qof(1) and f2o c (l) are expressed through the continuous fractions following from 
recursions: 

n ^ = (W O < 4.1V < 27 ) 

uj — [E F + 2L) jh — ujq — \lQ F {n + 1) 

2 

nil j 

n ^ = FT* al o 1 TTv ( 28 ) 

u — E c /n — u — Uo c (n + 1) 

= £f^0, ^al = £^0- (29) 

(2) In the one-phonon vibronic region 

ai = hLU- (E F + 2L) - [a 2 + MT] , (30) 

// . // . ^a^al / Q1 x 

«i2 = e c + e h + — — 0-3, (31) 

UJ 2 

a 2 = hj-E c -^[a 4 + M F T] , (32) 
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H[u-uj - fii c (l)] - E c 

a 2 = a, 

m 

a 3 = (1/m) «V + e£<7i) , 
1 



and 



0-4 



ft[u; -o; - fiic(l)] - -Ec 



1 + - 
m l 



«) 2 ct + (2< + 



T = — 



a 

m 



2 (e" N2 



(33) 
(34) 

(35) 
(36) 



m = 2 {L [ft [w - u - ^ic(l)] - E c ] + 2e'^} , (37) 
t = (1/m) I {ft [u - uj - fliF(l)] - £f} {ft [w - w - fiic(l)] - #c} 



- 2 



CO 2 + (^) 2 



D = 1 - Mpo~2 - Mcr 4 - M F MT, 

-l/v 7 ^ 7 ! if li + VtMkl, 



(7 



(38) 
(39) 

(40) 



l/v 7 ^! if t + <1, 



cti = to - 1. (41) 
The functions Qif(1) and Oi c (l) also represent continuous fractions from recursions: 



Qi c (n) 



nut, 



uj — Ep/h — (n + 1)ujq — ^if(^ + 1) 
null 



uj — E c /h — (n + 1)uj — Q\ c (n + 1) 
Finally, the functions M and M F can be expressed as follows: 
M/ft = fi 0c (2)-fi lc (l), 

M F /ft = n F(2) -fiiF(i). 



(42) 
(43) 

(44) 
(45) 
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Table 1. Excitonic and vibrational parameters (in cm of the MePTCDI and 
PTCDA crystals. 





Ep E c 




L 




2ip 


7 


MePTCDI 
PTCDA 


17992 17346 
18 860 18 300 


1400 
1400 


355 
330 


-380 -137 
-48 -436 


36.8° 
82° 


68.1° 
143° 



4. Simulations of the excitonic and vibronic spectra of MePTCDI and 
PTCDA crystals 

In this section, we calculate the absorption spectra of the two crystals finding the 
imaginary parts of the components of the linear optical susceptibility, ( |22l) and ( l23l , 
at (2p F /v) = 1 and supposing an imaginary part equal to i5/h of the frequency u. We 
put the excitonic and vibrational parameters for the studied crystals as they have been 
fitted in Refs. [5] and [10] and used in our previous papers [HI [12], [13], [Hj, see table 1: 
The data for angle 7 and ratio a = \pct/pf\ have been calculated in [6] using quantum 
chemical evaluations. We use the values a = 0.1 for the PTCDA crystal and a = 0.135 
for the MePTCDI. The data for angle 2(p are derived from the crystal structure (for 
MePTCDI see [21]). 

The aforementioned parameters are permanent in our calculations. We vary the 
values of the following parameters (intending to observe their impact on the absorption 
spectra and find a better similarity with the experimental absorption spectra [5], [6], [7]): 

(i) The excitonic damping quantity 5 which varies from 80 cm -1 to 500 cm -1 . 

(ii) The linear exciton-phonon coupling parameters £f and £ which vary near the 
values 

£ F = 0.82, f = for PTCDA (see 0), 

£ F = 0.88, £ = for MePTCDI (see [7]). 

In calculating continuous fractions ( 1271) . ( 1281) and ( ]42l) . (l4~3l) we take twenty steps. 

4.1. MePTCDI 

The general structure of the absorption spectra at 5 = 300 cm -1 can be seen in figured] 
in which the red curve is calculated using 'excitonic' formulas (I24 |) -(l26 |) (in the following 
denoted as 'exc' program) and the green curve by using one-phonon vibronic formulas 
(l30l)-( !32l) (the corresponding program is denoted as 'lp')- The approximate boundary 
between the two regions is 17800 cm" 1 . Two absorption maxima associated with the 
two types of excitons (FEs and CTEs) appear in the excitonic spectra below 17800 
cm -1 , as well as in the vibronic spectra (green curve). Our calculations are valid for the 
one-phonon vibronic spectrum, approximately between 17800 cm" 1 and 19100 cm -1 
but the higher vibronics also appear in the calculated spectra and they show splitting 
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associated with a FE-CTEs mixing (see the doublet near 19 500-20 000 cm -1 ). The 
shape of the experimental absorption curve, see figures 2 and 3 from Ref. [6], can be 
combined approximately using the excitonic curve (up to 17800 cm -1 ) and the green 
curve above that frequency. 

The impact of the FE-CTEs mixing can be seen more clearly by comparing the 
two curves in figure [2] calculated with the 'exc' program for two different values of the 
damping 5. The smaller value 5 = 80 cm -1 generates sharp absorption maxima in the 
excitonic and one-phonon vibronic regions. 

The one-phonon vibronic spectrum calculated with the 'lp' program is presented 
in figure [3j The maximum near 18 800 cm -1 at 5 = 300 cm -1 would be relatively flat 
but still resolvable. 

The impact of the CTEs transition dipole moment is illustrated in figure HI The two 
curves — the red calculated with the CTEs contribution, and the green one calculated 
without this contribution [a = in formulas (l2~TT) - (!23l ] — are relatively close to each 
other. The strongest impact can be measured near lower excitonic peak generated by 
the CTEs level E c = 17346 cm" 1 . 

Figure [5] demonstrates the changes of the absorption lineshape in the vibronic 
spectra depending on the values of the damping 5. The bigger value of S (the green 
curve) makes the absorption curve flat. The half- width of the highest red maximum 
approximately coincides with the magnitude of 5 (= 300 cm -1 ) which confirms its one- 
particle nature (bound exciton-phonon state). 

The two curve in figure [6] calculating by using the 'exc' program for 5 = 300 cm -1 
differ by the mutual positions of the two excitonic levels. The red curve corresponds 
to a lower position of the CTEs level E c < Ep, whereas the green curve corresponds 
to the opposite situation E F < E c (the magnitudes of E c and E F are exchanged). The 
green curve exhibits a strong domination of the lower maxima associated with the FE 
level. The two curves can be approximated with five Lorentz maxima in the spectral 
region of 15 000-22 000 cm -1 [6]. Comparing the lineshapes of these two curves with the 
experimental curve in Ref. [6], we cannot make a hypothesis which possibility is more 
probable. We prefer the fitting from Refs. j6] and [10] (E c < Ep), however, another 
choice of the mixing arrangements is also allowed. 

The importance of the FE-CTEs mixing can be seen in figure [7] in which the 
green curves represent a pure FE absorption (e" = e'^ = 0). The complexity of the 
experimental absorption curves with five Lorentz maxima in the studied spectral region, 
see figures 2 and 3 in Ref. [6] , cannot be understood on the basis of simple Frenkel-exciton 
model. Contrary, the model of mixed FE-CTEs and their vibronic satellites can be the 
basis of adequate simulations of the absorption curves (red curves). 

Figure [8] contains the calculated absorption curves at 5 = 300 cm" 1 (with the 'exc' 
program) for three different values of the linear exciton-phonon coupling constant £f- 
The values of £f which we put in calculations can be considered as a characteristic of 
the intermediate coupling, but the lineshape of the absorption curves is quite different. 
Obviously only the green curve exhibits five absorption maxima despite of the four 
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maxima and the two saddle parts of the blue curve. The value £p = 0.88 which we use 
in calculating other cases (figures [IHTj) seems to be the most probable. 

The model of the FE-CTEs mixing reproduces the general structure of the 
absorption spectra in the excitonic and vibronic regions of the MePTCDI crystal. Our 
calculation confirm the correct choice of the excitonic and vibrational parameters in the 
fitting procedure implemented in [61 [TO]. We establish as the most probable value of the 
exciton damping 5 = 300 cm -1 . For us, however, one open question still exists, notably 
the mutual positions of the excitonic levels of FE and CTEs lines. 

4.2. PTCDA 

The experimental linear absorption lines for the PTCDA crystal are broader than those 
of the MePTCDI (see [6]). Thus the structure of the PTCDA spectra does not exhibit 
many details due to the FE-CTEs mixing. 

In figure [9] the boundary between the excitonic and vibronic spectra is around 
18 800 cm -1 and the FE-CTEs splitting can be observed in the red 'exc' curve only. 
The splitting in the vibronic satellites is practically unresolvable. 

Figure [10] demonstrates relatively a weak impact of the FE-CTEs mixing on the 
vibronic spectra where the pure FE satellites (green colour) and the satellites of the 
mixed excitons (red colour) lie very close. 

The absorption curves in figure [TT1 have been calculated with £p = 1. The excitonic 
and vibronic lines at 5 = 300 cm -1 demonstrate doublet-like structure confirmed on the 
blue curve corresponding to a non-realistic low damping 5 = 80 cm -1 . The red and 
green curves are similar to the lineshape of the experimental curves. 

The real width of the excitonic curve is about 5 = 500 cm" 1 (figure [T2|) . The 
absorption above 18 500 cm -1 with £p = 0.82 (red curve) can be described with one 
maximum near 19 500 cm -1 and its higher vibronic replicas. The blue curve calculated 
with £p = 1.1 is more similar to the experimental curve (see [5]) whose lineshape covers 
several vibronic maxima. We accept as a very good simulation of the absorption spectra 
namely the case of 5 = 500 cm" 1 , £p = 1-1- It can be seen in figure [TBI that the excitonic 
and first two vibronic maxima represent doublets whose splitting results from the FE- 
CTEs mixing, but this doublet structure is hidden due to the big damping. The blue 
curve calculated at 5 = 100 cm -1 shows the true structure of the wide absorption 
maxima. 

Figure [TH shows insensitivity of the absorption lineshape in the case of a mutual 
replacement of the two excitonic levels E c < Ep and Ep < E c . Obviously the absorption 
curves are shifted, however, their lineshapes are very similar. 

So, in the PTCDA crystal the large excitonic linewidth which we estimate to be 
near 5 = 500 cm" 1 covers the effects of the FE-CTEs mixing. Another result of the 
simulations of the absorption spectra is the hypothesis for a higher value of the constant 
£f of the linear exciton-phonon coupling evaluated up to now as £f = 0.82 [SI Ej. Our 
calculations show that possible values of £p should be 1 or 1.1. 
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5. Conclusions 

Our model for the linear absorption spectra of the molecular crystals includes the 
following parameters: 

• excitonic levels Ep and E c , as well as the vibrational frequency u of the 
intramolecular mode, 

• parameters e" and of the FE-CTEs mixing, 

• constants £f and £ of the linear exciton-phonon coupling, 

• the width S of the excitonic linewidth, 

• angles 2ip, 7 and the ratio a of the CTEs and FE transition dipole moments. 

Practically all parameters have been introduced and fitted for the MePTCDI and 
PTCDA crystals in previous papers [5j El [7J [10] . In our study we apply the complex 
vibronic approach in calculating the linear optical susceptibility and its imaginary part 
which is a factor in the absorption coefficient in the excitonic and one-phonon vibronic 
spectra. Higher vibronics — with two, three phonons — also have been demonstrated 
in our calculations. The main goal of our model is to simulate the lineshape in the 
absorption region of 15 000-23 000 cm -1 of the aforementioned crystals, and to find out 
the adequate values of the excitonic linewidth, exciton-phonon coupling parameters and 
so on. The better coincidence between the calculated and measured absorption curves 
can be achieved by fitting all the parameters of the model, especially the positions and 
distances between the excitonic levels. The inclusion of external phonons can also cause 
some changes in the simulations [18], as well as the inclusion of several intramolecular 
vibrational modes. 

We stress again the main conclusions which concern the two crystals: 

(i) The excitonic linewidth of the MePTCDI crystal, according to our simulations, is 
approximately 5 ~ 300 cm -1 . The FE-phonon linear coupling coefficient £p has 
been estimated correctly in the previous papers [HE] as £p ~ 0.88. The absorption 
spectra depend strongly on the mutual position of the two excitonic levels Ep and 
E c , but the simulations of the absorption spectra in our paper do not give the 
opportunity to choose one of the two possibilities: Ep < E c or E c < Ep. 

(ii) The excitonic linewidth of the PTCDA crystal can be evaluated as 5 ~ 500 cm -1 . 
That is why the effect of the FE-CTEs mixing, being covered by the wide absorption 
maxima, are weakly expressed than in the MePTCDI crystal. However, a very 
probable conclusion from our calculations may be the stronger linear exciton- 
phonon coupling (£p ~ 1 or 1.1 instead of 0.82). 

A supposition for FE-CTEs mixing is necessary for the interpretation of the 
excitonic and vibronic spectra of the MePTCDI crystal. Concerning the PTCDA crystal 
our calculations suggest two possible models: (a) the FE-CTEs model studied in this 
paper, and (b) the model of pure FE and its vibronic satellites assuming wide excitonic 
levels. 
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Our model can be applied in the interpretation of other one-component molecular 
stacks (crystals) and it can be more effective in the systems with narrower excitonic 
absorption lines. 
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Figure 1. Linear absorption of the McPTCDI crystal at S = 300 cm -1 and £f = 0.88. 
Red curve is calculated for the excitonic region [using formulas (|24[> — (f^^|) ] and the green 
curve for one-phonon vibronic spectra [formulas (|30[) (1321) ] . 
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Figure 2. Absorption spectra calculated with program 'exc' for the following values 
of the excitonic damping: red curve with S = 80 cm -1 and the green one with 8 = 350 
cm -1 . 
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Figure 3. Onc-phonon vibronic absorption spectra at £f = 0.88, S — 80 cm" 1 (red) 
and S = 350 cm" 1 (green). 
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Figure 4. Impact of the CTEs-transition dipolc moment expressed through the ratio 
a. Red curve for a = 0.135 and the green one for a = 0. 
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Figure 5. Vibronic absorption at various values of S: red curve for 8 = 300 cm -1 and 
the green curve for S = 400 cm" 1 . 
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Figure 6. Comparison of absorption curves calculated with data from table 1 [red 
and green curves calculated for E-p = 17 346 cm" 1 and E c = 17 992 cm" 1 , respectively 
(interchange of the excitonic levels)]. 
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Figure 7. Linear absorption spectra with FE-CTEs mixing (red curves) and pure 
Frenkel exciton spectra (green curves at e" = = 0), 6 = 300 cm -1 : (a) absorption 
spectra calculated with 'exc' and (b) 'lp' absorption curves. 




Figure 8. Linear absorption spectra 'exc' at S = 300 cm" 1 calculated for the following 
values of the FE-phonon coupling constant £f : re d curve for £p = 0.7, green curve for 
£f = 0.88, and blue curve for £p = 1. 
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Figure 9. Linear absorption of the PTCDA crystal at £f = 0.82 and 6 = 400 cm 1 . 
Red curve for the 'exc' calculations and the green one for the 'lp' calculations. 
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Figure 10. Linear absorption of PTCDA ('lp' program) at £f = 0.82 and 5 = 400 
cm -1 with FE-CTEs mixing (red curve) and without it (green curve with e" = e[[ = 0). 
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Figure 11. Linear absorption of PTCDA at £f = 1- Red curve ('exc') at 5 = 300 
cm -1 , green curve ('lp') at 6 — 300 cm -1 , and blue curve ('lp') at 5 = 80 cm" 1 . 




Figure 12. Linear absorption of PTCDA at S = 500 cm 1 . Red curve ('lp') 
corresponds to £f = 0.82 and the blue curve ('lp') — to £p = 1-1. 
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Figure 13. Linear absorption of PTCDA at £f = 1-1. Red curve ('exc') and green 
curve ('lp') correspond to 8 = 500 cm -1 , while the blue curve ('exc') — to S = 100 
cm -1 . 



0.001 




- 



Figure 14. Linear absorption of PTCDA at £f = LI and S = 500 cm -1 ('lp' 
program). The two excitonic levels are interchanged and the red curve corresponds to 
the values in table 1 (E c < Ep) whereas for the green curve Ep < E c . 
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